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RATE  COMPLEMENTARY  ENERGY  PRINCIPLES; 

FINITE  STRAIN  PLASTICITY  PROBLEMS;  AND  FINITE  ELEMENTS 


Satya  N.  Atluri 

School  of  Engineering  Science  and  Mechanics  ^ 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332,  USA 


ABSTRACT: 

Complementary  energy  theorems,  for  the  rate  problem  of  finite  strain  clas- 
sical elasto-plasticity,  in  both  Updated  and  Total  Lagrangean  rate  forms  with 
alternate  stress-rates  and  conjugate  measures  of  strain  rate,  are  studied  from 
the  point  of  view  of  their  application  in  finite  element  schemes.  Two  new 
complementary  theorems,  in  Updated  and  Total  Lagrangean  forms  respectively, 
are  proposed.  The  relative  merits  of  these  in  application  to  finite-strain 
elasto-plastic  stress  analysis,  in  the  treatment  of  near- incompressibility  at 
large  plastic  flow,  and  in  treatment  of  buckling  problems,  are  briefly  discussed. 


1.  COMPLEMENTARY  PRINCIPLES  IN  UPDATED  IAGRANGEAN  (UL)  FORM:  RATE  PROBLEM  OF 
FINITE  STRAIN  CLASSICAL  ELASTO-PIASTICITY : 

In  this  case,  the  current  state  C^,  wherein  the  carters ian  spatial  co- 

N N 

ordinates  of  a particle  are  y and  the  true  (Cauchy)  stress  is  £ , is  used  as  a 
reference  state  for  measuring  the  rates  of  variables,  in  going  from  CN  to  CN+^ 

If  d is  the  corotational  rate  of  Kirchhoff  stress  a (=  J t ; J is  the  deter- 
minant  of  the  Jacobian);  £ is  the  UL  strain  rate  {2  1/2  [(VNu)  + (V^u)^]] ; £ the 
UL  spin  rate  {s  1/2  (V^u)T  - (V^ju)]};  the  gradient  operator  in  C^  (=  .e^  d(  )/ 
dyt);  u is  the  rate  of  displacement;  £ and  £,  respectively,  are  the  rates  of 
second  and  first  Piola-Kirchhof f stresses  referred  to  and  measured  per  unit  area 


of  C„;  it  is  well-known  that: 
N 


. . N N . • .*  . N N. 

s - u - e.T  - t . e ; t = d - d.T  - t .(I) 

**  r>j  r*j  r>J  rsj  m 


(1.1,2) 


Further,  from  the  relation  for  the  Jaumann  stress  tensor  r [l],  viz.,  r = 1/2 

/v  M 

T X 

[t  . a + a . £ ] ; we  obtain  the  equation  for  the  UL  rate  of  r as : 

f = 1/2 (t  + tN  . (I)  + * . tn  + tT)  = d*  - 1/2 (d  . tN  + tN  . i)  (1.3a,b) 
As  noted  by  Hill  [2],  the  constitutive  equation  of  a bilinear  form,  for  work- 


hardening elastic-plastic  materials,  in  terms  of  d and  2 can  be  written  as: 


2* ' Wi/w ' I <Wu>: 


(1.4a, b) 


Where  L^^,  O',  g,  and  X^  are  as  defined  in  [2].  From  Eqs . (1. 1-1.4)  it  is  seen 


that , 


s = dW/de  ; W = V - tN:  (e  . e) 

/V  A/'  A*  >V  /V 


t = 90/*/;  U = V - tN:  (£.£)  + (1/2)tN:  (&T  . e) 


(1. 5a,b) 


(1.6a,b) 


f = 9Q/de  ; Q = V - (1/2)T  : (£  . £)  (1.7a,b) 

N X 

where  e = (V  u)  . It  can  then  be  shown  [3,4]  that  the  UL  rate  form  of  a com- 
plementary energy  theorem,  with  £ and  u as  variables,  can  be  stated  as  the 


I 


stationarity  condition  of  the  functional. 


*^(s;u)  = r { -S*(s)  - (1/2)tN:  [ Q7Nu)  . (VNu)T]}dv  + J*  t . uds  (1.8) 

c ~ w ~ ~ ~ <?-” 

N “un 
with  the  a priori  constraints: 

Linear  Momentum  Balance  (LMB) : VN  . [s  + tN  . (VNu)]  + pNB  =0  (1,9) 

• T 

Angular  Momentum  Balance  (AMB)  : £ = s,  (1.10) 

Traction  Boundary  Condition  (TBC)  : n . [j^  + £ • (V^u)]  = _t  = _t  at  S (1.11) 

N 

and  S is  pressumed  to  be  obtained  through  the  contact  trans format ion , 

s**<i>  = i’i  • w(£)  d.i2) 

and  with  the  Euler-Lagrange  Equations  (ELE)  and  Natural  Boundary  Conditions  (NBC): 
(i)  compatibility,  £ = 1/2  [(V^u)  + (V^u)'*']  and  (ii)  Displacement  B.C  (DBC) , 
u = u at  S . The  constraint,  of  TBC  is,  in  general,  impossible  to  satisfy  a 

“n 

priori;  whereas,  it  may,  in  some  special  circumstances,  be  possible  to  satisfy 

the  constraints  (1.9  and  1.10).  For  instance,  first  noting  that  £ + £**  • (£Nu) 

* ’ N ' N’  ' N N 

= t,  one  may  assume  a ^t  that  satisfies  UiB  V . t + p B = 0 (by  setting  Jt  = V xY(y  ) 

+ t?  where  ^ is  any  particular  solution  such  that  ~ P**B)  and  also  assume  an 

arbitrary  but  symmetric  s field  that  hence  obeys  the  AMB  condition.  One  can  then 

eliminate  displacements  from  Eq.  (1.8)  by  setting 


VNU  = (tV1  . [VN  x Y + tP  - s] 


(1.13) 


Thus,  provided  that  the  principal  values  of  £ are  not  zero  (which  is  a rather 
special  circumstance)  one  can  express  Eq.  (1.8)  in  terms  of  stress-rates  only. 

On  the  other  hand  one  can  start  out  assuming  u and  then  set, 

s = curl  curl  A + s^;  where  . s*5  = -p^B  - . [t^  . (V^u)J  (1.14) 

where  A are  the  well-known  Maxwell-Morera-Beltrami  stress  functions  and 

rss  ^ 

is  any  symmetric  particular  solution,  such  as,  for  instance,  given  by: 

s*\  = 0 (i*j);  = Jy  ^[-pNBi  “ d?i  <no  sum  on  l>  (1,15) 


A 


However  the  question  of  completeness  of  ^sP  (or  the  consquences  of  a priori  set- 
ting sP.  =0,  i 4 j)  remains  to  be  answered.  Moreover,  assuming  both  I and  u 
° 1 j ~ 

as  above  has  the  less  desirable  consequence  that  the  displacements  derived  from 
strain-integrals  corresponding  to  the  assumed  stresses  do  not,  in  general, 

agree  with  the  assumed  displacements.  Lastly,  one  can  first  assume  u and 
then  let. 


s = curl  curl  A - TN  . (VNu)  + sP;  VN  . sP  = -pNB 


(1.16) 


However,  the  chosen^  above  ceases  to  be  symmetric  and  the  AMB  is  violated. 
Eventhough  the  representations  as  in  Eqs . (1.14,15  & 16)  appear  to  be  less 
than  desirable  in  general,  use  of  such  have  been  made  in  structural  mechanics 
of  curved  beams  and  shells,  and  their  stability,  in  [5,6]. 

On  the  other  hand,  an  alternate  general  variational  theorem,  with ^t,  £, 
and  u as  variables,  can  be  stated  [3,4]  as  the  condition  of  stationarity  of 
the  functional: 

= J £u(e)  - pNjj  . u+tT:  [ (VNu)T  -J]}dv  - J t.uds  - J £ . (u-u) 


V, 


N 


which  leads  to  the  following  ELE  and  NBC: 


N 


U 


N 

(1.17) 


(LMB):  VN  . t + pNB  = 0;  (CC):  & * (VNu)T;  (DBC) : u = u at  S ((1.18,19 

“N 

(RCL):  2itJ  . 3,0/24^  - aijkl  - <*ta  + 


N . . . . 

" Tix(ejej  " fejx* 


(1.21) 


N.,T  N 


(AMB):  (V u)  . T + t = SYMM;  (TBC) : t = (n  . ^t)  - t at  S 


(1.22  & 


ds 


& 20) 


i 


23) 


1 1 


< 


■ 


The  AMB,  Eq.  (1.22),  is,  however,  embedded  in  the  structure  of  U when  U is 

defined  according  to  Eq.  (1.6b).  It  is  also  noted  that  the  RCL,  Eq.  (1.21), 

cannot,  in  general,  be  inverted  analytically.  Assuming  that  Eq.  (1.24)  can  be 

inverted  numerically,  one  can  achieve  a contact  transformation,  ^t  : £ - U 

. *.  • 

(e)  = T ^t) . If,  further,  one  satisfies  a priori  the  Eqs.  (1.18  & 23),  one 
obtains  a complementary  energy  principle  governed  by  the  stationarity  of: 


%2c('t)  = J*  “T * (t)  dv  + r t . u ds  (1.24) 

VM  S 

N u 

N 

which  has,  as  its  ELE  and  NBC,  Eqs.  (1.19  & 20).  This  principle  was  first 
stated  by  Hill  [7].  However,  the  AMB  condition,  Eq.  (1.22)  is  not  an  ELE.  Thus, 
for  the  principle  to  be  valid,  the  AMB  condition,  if  at  all  it  is  met,  must 

* • * . . * . T 

be  embedded  in  the  structure  of  T' ^t) ; ie.,  T* ^t)  must  be  such  that,  t + 0 T*/dt  ) 
• Z,  =~  + X * ©T’/d^t).  It  can  be  shown  [4]  that  this  is  not  true  in  general; 
and  thus  the  above  is  not  a valid  principle,  in  general. 


To  avoid  the  above  difficulties,  we  transform  the  functional  in  Eq.  (1.17) 
to  one  with  r;  u;  u>;  and  e as  variables.  First  we  note  that  e = e + ui ; and 
u = Q + (1/2)tN:  (ujT  . u>)  + tN:  <juT  . £)  . When  these  substitutions  are 
made,  we  rewrite  the  functional  in  Eq.  (1.17)  as: 

i>  & ■ J CO  + <1/2>iN:  • <s>  + xN;  <“>T  • V “ pN^-^  + 

VN 

+ tT:  [ (VNu)T  -£  - up 3 dv  - J t.  uds  - J t . (u  - u)ds  (1.25) 

~ ~ ~ S„  S 

°N  UN 

With  the  constraint  that  u>  must  be  skew-symmetric.  The  stationarity  of  the 
above  functional  leads  to  the  ELE  and  NBC: 


(RCL)-  SQAc  = (1/2)  C£  + f . w + <iT.  rN  + tT)  = i 

(LMB)—  VN  . t + pNB  = 0;  (AMB)-  tN  . is  + c . tN  + t = SYMM. 

(CC)-  (VNu)T  = e + (i;  (TBC)—  n"  . t = t at  S ; u = u at  S 
~ ~ ~ — a — — u 

N N 


(1.26) 

(1.27  & 28) 
(1.29,30  & 31) 


5 


It  is  noted  that  (AMB)  is  a clear  ELE  (corresponding  to  variation  inw)  and 
need  not  be  embedded  in  the  structure  of  Q.  If  we  establish  the  contact  trans- 
formation, that^r:  £ - Q = R ^r)  and  further  if  1MB  and  TBC,  Eqs.  (1.27  & 30), 
are  satisfied  a priori,  we  can  eliminate  u and  £ as  variables  from  Eq.  (1*25) 
and  obtain  the  following  functional  whose  stationarity  condition  leads  to  a 
complementary  energy  theorem: 


S2C<&&  = I C-  + (1/2)tk:  (&l  . <fi)  - tT.«:}dv+Js  (n*^)  . uds  (1.32) 

VN  • " " “N 

* • *T  N . .T  N 

Where,  by  definition,  ^r  = l/2^t  + ^t  + £ . to  + uj  . £ ) . In  the  above  functional 

the  constraints  are  the  LMB,  Eq.  (1.27)  (which  is  easy  to  meet  by  setting 
t = VN  x Y + tP;  tP  is  any  particular  solution  such  that  V*1  . tP  * -p**B),  and 
that  oo  is  a skew- symmetric  field  [which  is  easy  to  meet,  by  setting  u> 

The  ELE  and  NBC  are  the  AMB,  Eq.  (1.28);  CC,  Eq.  (1.29);  and  the  DBC,  Eq.  (1.31). 

In  as  much  the  constraints  in  the  principle  are  trivial  to  be  met  a priori,  and 
the  AMB  is  clearly  an  ELE,  the  above  complementary  principle  appears  to 


N 


,,T 


“ij 


eijkV] 


be  the  most  rational  and  consistent  for  practical  applications.  In  applying 

the  above  principle  in  conjunction  with  a finite  element  method,  the  initerelement 

* . + * • - 

traction  reciprocity  constraint,  (n  . t)  + (n  . t)  = 0 at  p ; can  be 

”n 

relaxed  a priori,  and  enforced  through  Lagrange  Multipliers.  In  doing  this, 
either  the  well-known  hybrid-stress  FEM  of  Pian  or  the  equilibrium  model  of  F. 
de  Veubeke  can  be  used. 


*It  is  worth  noting  that  the  above  principle  is  valid  for  both  isotropic  as  well 
as  anisotropic  materials.  In  the  case  of  isotropy,  it  can  be  shown  [4]  that  the 
AMB  reduces  to  the  constraint:  . d>  + t is  symmetric.  However,  even  for  iso- 

tropy, it  is  more  convenient  to  retain  thTs  AMB  condition  as  an  ELE  corresponding 
to  variations  in  d>.  Thus,  when  £ (that  satisfies  LMB),  and  a skew- symmetric  <I> 
are  assumed  (which,  however,  need  not  obey  the  AMB  even  for  isotropic  materiaTs), 
the  above  definition  for  r should  be  used  to  recover  the  AMB  as  an  Eli  for  either 
isotropic  or  anisotropic  materials. 


6 


! 


2.  TL  RATE  COMPLEMENTARY  PRINCIPLES:  RATE  PROBLEM  OF  FINITE  DEFORMATION 
ELASTO -PLASTICITY: 

In  this  case,  the  initial  (undeformed)  state  CQ  is  used  as  a reference 

N N 

state  in  all  subsequent  considerations.  Let  t and  s be  the  first  and 

second  Piola-Kirchhoff  stress  tensors  in  C„  as  referred  to  and  measured 

N 

per  unit  area  in  Cq;  and  let  ^t ' and ^s'  be  the  corresponding  TL  rates;  and 

let  E'  be  the  TL  rate  of  strain.  It  can  be  shown  that, 

**  • 

E'  = l/2(e’  + e,T  + e,T  . eN  + eNT  . e');[  where  e'  = (V°u)T;  eN  = (V°uN)T] 

“ <£N)T  * £ • [where  FN  » (V°£N)T]  (2.1a,b) 


N N.  -1  . .-N.-T,  , tN._,N.-1*  , JJT  . N ,T 

J(F)  . « . (F  ) ;t'=J(F)  t = s . F +s  ,e* 


(2.2  & 3) 


l/2[ tN  . a'  + a'T  . tm  + t'  . oN  + am  . t,T] 
********  **  ******  **  J 


(2.4) 


N N N 

Further  we  consider  the  polar  decompositions,  F = a . (I  4*  h ),  where  -f  h 

N N 

is  called  the  stretch  tensor;  and  e'  = a 1 . (I  + h ) + ty  . h' . The  rigid 

N NT  N 

rotation  tensors  satisfy  the  orthogonality  conditions,  a .a  = ^1;  and  a 
T NT 

a'  + a'  . a =0.  Because  of  Eqs.  (2.1b  & 2)  it  is  seen  that  ifW  is  a 
****** 

potential  forje,  a potential  W'  can  be  derived  for  ^ such  that,  ' = dW/dE'. 
Further,  because  of  Eqs.  (2.3  & 4)  it  is  seen  that: 


t'  = &U  7de'T;  U ' = w'  + (1/2)  sN:  (e,T  . e') 
«.■  **  9 ' **  **  ** 


(2.5a,b) 


r‘  = *Q7dh*;Q'  = W<£’)  + d/2)/:  Q}'  . h') 


(2.6a,b) 


From  the  above,  we  obtain  the  contact  transformations,  that, 

s’:  E’  - W’  = S’(s');  t,T:  e'  - U’  = T’(£’);  £*:  h’  - Q*  = R'(r')  (2.7a,b,c) 
Using  procedures  analogous  to  these  discussed  earlier  in  connection  with  the 
U1  rates  case,  and  as  elaborated  in  [4],  the  TL  rate  complementary  principle 
with  s' , and  u as  variables,  can  be  stated  as  the  stationarity  condition  of: 


7 


L 


^ = J fS'(V)  + (1/2)sN:  (je'T  . e'))dv  - J*  t'  . uds 


(2.8) 


TO*  ^ rp  © 

where,  <e'  = (V  u) ; and  with  the  constraints:  (AMB,  w'  = <s'  ; (IMB) : V 


fsN  . e,T  + s'  . FNT}  + p°B'  = 0;  (TBC) : n . {sN  . e ,T  + s ' . f”1}  - t'  = t'  at 

44  44  44  ^ r 44  0+4  0+4  44  J — 


. The  ELE  and  NBC  resulting  from  the  stationarity  of  the  above  functional 
o 


are:  CC,  Eq.  (2.1a);  and  DBC,  u = u at  Sy  . The  possible  ways  to  satisfy  LMB 


a priori,  are:  (i)  to  choose  ^t ' that  satisifes  LMB  and  in  addition  to  choose  a 

NT  NT 

symmetric  s',  thereby  elminate  e 1 as  a variable,  as:  s . e'  + s ' . f”  = t* 

0*4  *4  *4  *4  *4  *4  *4 


= v”  x Y + t'P;  (Vu  . t"*  = -p"B')  and  e 

— 44  44  7 X 44  r — • ' 44 


P - NnI, 


,T  . N. -1 


r1  . [V°  x Y + t'P  - s'  . FNT] 

* — 44  44  44  44 


which  requires  that  the  principal  values  of  £ be  non  zero.  On  the  other  hand, 
one  can  choose  u in  addition  to  t'  that  satisfies  1MB,  and  derived  s'  from 

— W 44 

them  as:  s'  = [v°  x Y + t'P  - sN  . e,T]  . (FNT) . ^ However,  s'  is  then  unsymmetric 
in  general  and  hence  violates  (AMB).  In  spite  of  this,  there  may  be  situations 
in  structural  mechanics  where  the  above  representation  may  be  satisfactory; 
however,  with  the  same  drawback,  as  discussed  earlier,  in  assuming  both  stresses 
and  displacements. 

It  can  be  shown  [6],  as  in  the  UL  rate  case,  the  TL  complementary  principle, 
analogous  to  Eq.  (1.24),  in  terms  of  ^t'  alone,  is  invalid  in  general.  Finally, 
the  TL  rate  complementary  principle  for  finite  strain  elasto-plasticity,  with 
t'  and  a'  as  variables,  which  is  the  counterpart  of  that  stated  in  Eq.  (1.32), 
can  be  stated  as  the  stationarity  of  the  functional: 


«&■*•>  -JC-  R*'<r')-fT*  [£’  . a + hN)]  - 1/2 /T:  & . g”  • 

V 

o 

a'  • tt  + h,N)]  } dv  + J t’  .uds  (2.9) 


Xi/here  ^r'  is  defined  as  in  Eq.  (2.4),  and  with  the  constraints:  (i)  V°  . t'  + 

p °B ' = 0,  which  is  easy  to  meet  by  setting  ^t'  = V°  x Y + ^t'P;  (ii)  crN  . a,T  + 
NT 

a’  . a =0  which  is  also  easy  to  satisfy  by  writing  = (*.^(9^)  where  9^ 


are  Euler-angles  of  rigid  rotation:  and  (iii)  n . t'  = t at  S . The  ELE 

- ~ — o 

N N N i-  N N ° 

and  NBC  are  (i)  AMB-*  h1  . t .a  + (h  + I)  . (t1  . a + t . or')  is  symmetric; 
(ii)  CO  (V°u)T  = a 1 . (^1  + h^)  + aN  . h'  and  the  TBC.  As  before,  the  above 
TL  rate  complementary  principle  is  believed  to  be  the  most  consistent  for 
practical  application. 

The  consistent  complementary  principles,  as  stated  in  Eqs.  (1.32  and  2.9) 
are  useful  in  (i)  obtaining  better  stress  solutions  in  a numerical  analysis  of 
finite  strain  elasto-plasticity  and  (ii)  in  treating  situations  of  near  in- 
compressibility as  many  arise  at  large  plastic  flow  and  (iii)  in  analysis  of 
stability  problems.  Numerical  evidence  for  (i)  and  (ii)  appears  to  provided 
in  the  finite  elasticity  solutions  for  compressible  as  well  as  incompressible 
materials  reported  in  [8,9].  In  the  stability  problems,  it  is  seen  that  the 
functional , 


h 


I 

i 

•! 


"I  {*'(£')  fe'  • « +£S)]  + 0/2)/1:  [«'  .&m 


(I  + hN)]]dv 


<V  /V 


(2.10) 


with  the  constraints  V°  . t'  = 0,  aN  . ,T  + a'  . <*NT  = 0;  and  n . t'  = 0 at  S 

— rv  " /s*  ew  rv  /v  7 — ^ q 

o 

attains  a zero  stationary  value  at  neutral  equilibrium.  In  an  FEM  procedure 
this  translates  into  the  criterion  that  the  current  stiffness  matrix  becomes  singular. 
If  the  prebuckling  state,  say  Co>  is  linear,  one  can  make  the  usual  approxima- 
tions: h = 0;  or  = I;  and  r'  = 1/2(1 ' + t ,T  + cr ,T  . t°  + t°  . ot ' ) : and  Eq . (2.10) 
becomes , 


. \ 


2 

= .f  ^'CS')  +£'  : 2'  + (1/2HoT:  (»'  .£')}dv  (2.11) 

o 

with  the  constraints  a'  = <v,T;  v°  . t'  = 0;  n . t'  = 0 at  , and  thus  leads 

, o 

to  an  eigm  value  problem. 
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